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We report the creation and real-space observation of magnetic structures with well-defined topo- 
logical properties and a lateral size as low as about 150 nm. They are generated in a thin ferri- 
^■^ I magnetic film by ultrashort single optical laser pulses. Thanks to their topological properties, such 

r^ , structures can be classified as Skyrmions of a particular type that does not require an externally ap- 

Q ■ plied magnetic field for stabilization. Besides Skyrmions, we are able to generate magnetic features 

O ' with topological characteristics that can be tuned by changing the laser fluence. The stability of 

such features is accounted for by an analytical model based on the interplay between the exchange 
and the magnetic dipole-dipole interactions. 

>■ 

]^ \ Skyrmions [1] are particle-like solutions of wave equations characterized by a topological index which is conserved 

in time and plays the important role of a quantum number for particle states in the corresponding field theory. One 
of the most interesting characteristics of such topological states of matter resides in the robustness that they hold 
with respect to perturbations and disorder. 

In magnetic materials Skyrmions emerge as soliton-like excitations that cannot be traced back to the ground 
ferromagnetic state by continuous deformations of the local magnetization field. So far they have been observed in 
materials with perpendicular magnetic anisotropy in connection with one of the following stabilization mechanisms: 
(i) four-spin exchange, leading to the formation of a Skyrmion lattice with each Skyrmion extending over a few lattice 
sites [2|; (ii) the Dzyaloshinskii-Moriya interaction (DMI) [3, 4], which is active in noncentrosymmetric helimagnets, 
where Skyrmions have typical dimensions of few tens of nanometers [5-8]; (iii) the long-range dipole-dipole interaction 
;h ■ (DDI) [9] stabilizing Skyrmions with a typical lateral size of the order of 1 fiva [10-12] . It has been recently found 
^_' that the DDI can induce a larger and more complex variety of magnetic textures with respect to the ones that are 
observed in association with four-spin exchange and the DMI [9|. 

Despite the native state of thin films of magnetic materials with perpendicular anisotropy is characterized by stripe 
domain patterns, Skyrmions can be generated when a sufficiently high magnetic field is applied perpendicular to the 
film surface. In this case, the stripes undergo a transition to a hexagonal Skyrmion lattice [THIQI, [12|. Skyrmions, 
however, can also be produced individually, as demonstrated by a laser-induced magnetization reversal experiment on 
a TbFeCo thin film [13], which still required a constant magnetic field to avoid the collapse of the domains produced 
upon illumination by single laser pulses. It has been proposed LM] that such magnetic features could be stabilized 
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FIG. 1: (Color online) (a, b) Near- field Faraday rotation map showing magnetic domains induced in a thin TbFeCo film after 
single laser pulse irradiation (energy density ^ 5 mJ/cm^) in film areas showing opposite out-of plane magnetizations, (c, d) 
Corresponding Skyrmion spin textures acccording to Eq. Q (with Oo = 7r/2). (e, f) Profile (dots) of the Faraday rotation 
measured along the lines in (b) and (c), respectively. The lines in (e) and (f) are fits performed with the trial function describing 
the Skyrmions discussed in the text [Ri = 200 nm for the structure in (a) and Ri — 235 nm for the structure in (b)]. 



in the form of small cylindrical domains thanks to the frustration among the exchange interaction, the magnetic 
dipole-dipole coupling, and the external magnetic field. 

We report the creation of magnetic Skyrmions in a ferrimagnetic thin film generated without any external field 
with the help of an ultrashort laser pulse. In addition, by varying the laser fiuence we are able to control the size 
of the Skyrmions and their topological quantum number, generating either Skyrmions or even states of a Skyrmion- 
antiSkyrmion pair, which we named Skyrmionium in analogy to similar bound particle/anti-particle states. In the 
following we demonstrate that a model including exchange coupling, single-ion magneto-crystal anisotropy and mag- 
netic dipole-dipole interactions can indeed have metastable solutions with a nonvanishing topological index. We 
show that the stabilizing action of the externally applied reversed magnetic field exploited in Refs. |l3| andH can be 
effectively replaced by the dipolar field from the magnetic film itself. 

The sample is a 20 nm-thick amorphous thin alloy film of Tb22Fe69Co9 which was deposited by magnetron sputtering 
on a Si3N4(5 nm)/AlTi(10 nm)/glass substrate and then capped with a 60 nm-thick layer of Si3N4. The thin metallic 
layer of AlTi acts as a heat sink. The film exhibits out-of-plane anisotropy with a coercive field of 1 T and a Faraday 
rotation between the two opposite magnetizations as high as 0.3° at A=632.8 nm. 

The sample was prepared by irradiation, at room temperature, with 150 fs circularly polarized single laser pulses 
(spot diameter ^ 2 /im, A = 800 nm), to induce ultrafast magnetization reversal [i5|. Permanent sample damage 
occurs for laser fiuences above 15 mJ/cm^. At variance with Ref. 13 and the works reporting DMI-stabilized Skyrmions 
[5-8], no external magnetic field was applied during this step. Finally, we mapped the out-of-plane component of the 
local magnetization with sub-diffraction resolution by measuring [16j the Faraday rotation of linearly polarized light 
(A = 635 nm) transmitted through the film from a polarization-conserving hollow pyramid tip (aperture diameter 
?^ 100 nm) [17] that was scanned over the sample. 

As shown in Fig. [TJ magnetic domains are obtained after single laser pulse irradiation with a pulse energy density 
of about 5 mJ/cm^. Within experimental accuracy, these sub-/im magnetic domains display a local magnetization in 
their center that is completely reversed with respect to the surrounding material. Their lifetime is very long as they 
have not shown any noticeable change more than 12 months after the first measurements. 

The mechanism stabilizing these structures resides in the interplay among the exchange coupling J^x, the unidirec- 
tional single- ion anisotropy J^a, and the magnetic dipole-dipole interaction J^d contributions to the total Hamiltonian 
H = Hx + Ha + i^D + Hz^ with Hz being the Zeeman energy (see the Supplementary Material). The value of the 
uniaxial anisotropy term Ha needs to be chosen in order to also account for the presence of domain wall pinning 
centers, which are responsible for the high coercive fields characterizing relatively hard magnetic materials such as 
TbFeCo. Both magneto-crystal anisotropy and pinning centers, in fact, have a similar effect in terms of creating 
an energy barrier for magnetization reversal. This approach consisting in treating pinning centers as a source of 
anisotropy is valid as long as i) pinning sites are uniformly distributed and ii) a single Skyrmion comprises many of 
them. It is reasonable to assume that these two conditions apply to our case since the TbFeCo film is amorphous and 
the solitons that we observe have dimensions of the order of hundreds of nanometers. 

In order to evaluate the existence of stable solutions to the integro- differential equations that stem from the Hamilto- 



nian iJ, suppose there exists a soliton solution n°(r). We can calculate the contributions to the total energy associated 
with n^(r) and denote them as E^ (exchange), E^ (anisotropy), E^ (DDI) and E^ (Zeeman). Now consider a new 
configuration, n(r) = n^(7^r). By evaluating the scaling properties of Hx, i^A, ^D, and Hz, we find the energy of the 
new configuration to be 

EUV) =Ei + r^El - r'\El\ + V-'El (1) 

This function has a unique minimum at 

Therefore the problem is not scale-invariant and, according to the Derrick-Hobart theorem [l9|,|20|, may have a unique 
soliton solution characterized by 77 = 1 or 2 (£^^ -\- E^^ = I^dI- ^^ ^^ apparent that the presence of perpendicular 
anisotropy can stabilize a soliton even without the Zeeman term, as already suggested by Clarke et al. [21|. Although 
the role of the DDI has been questioned [22|, Eq. (J2j) demonstrates that the DDI is an essential requisite for the 
stabilization of the soliton. 

We proceed by examining cylindrically symmetric general solutions of the nonlinear sigma model derived from the 
expression of the Hamiltonian H having n^ = cr(^), with r being the distance from the soliton center and (j(r) a 
trial function, —1 < a{r) < +1. The Pontryagin number Q [18] of this spin structure is readily calculated to be 
Q = [cr(oo) — cr(0)]/2 = +1, while the total energy is finite and positive. These two conditions together with the 
Derrick-Hobart theorem allow us to assert the existence of a function that minimizes the total energy, even without 
knowing its explicit form. The winding of the spins in the structure is described by a zero-energy mode parameter ^0 
(see Supplementary Material), which affects the total energy only through small corrections to the DDI term and has 
a negligible infiuence on the soliton size. 

Since it is not possible to find analytical solutions, we just consider a set of reasonable trial functions and determine 
the one that minimizes the total energy. The line profiles in Fig. [1] reporting the spatial variation of the Faraday 
rotation indicate that a{r) can be chosen as a Gaussian function, 

a(r) = l-2e-^'/^i, (3) 

resulting in a spin texture with Q = 1 (in the following we will adopt the notation Rq to indicate the size parameter 
for solitons with Pontryagin number equal to Q). 

The parameters Ri and ^0 can be determined by minimizing the total energy Hx + Ha + H^ + Hz of the soliton 
(see the Supplementary Material), giving 

i^6m 

' nT/e-2n/df^2^Az/a^' ^ ' 

df being the thickness of the film, F the exchange constant, ^ the single- ion uniaxial anisotropy, Q the DDI strength, 
a the lattice constant, and Az the Zeeman energy (see the Supplementary Material). The ground state has spins 
forming a Bloch-type magnetic domain wall (^0 = i7r/2) separating the soliton core from the surrounding medium. 
This soliton has a finite radius and is stable even without an external field (Az = 0) in a sample characterized by 
finite uniaxial anisotropy (^ 7^ 00). Our analysis of the stability of spin arrangements consistent with Fig. [1] and their 
nonvanishing Pontryagin number allow us to consider them as genuine Skyrmions. 

Another feature differentiating such Skyrmions from the disc-type ones discussed in Ref. 14 is that the latter are 
stable only above a minimum size, while the size of the former can be shrunk without limit by increasing the magnetic 
field. By using typical sample parameters [14, 23], namely F = 3.5 x 10^^ J, ^ = 25 nm, and a = 0.3 nm, we reproduce 
the experimental value of the Skyrmion size in zero field, which is about 200 nm (see Fig. [1]) for (^ = 1.6 x 10^^ J/m, 
which would correspond to a saturation magnetization of about 56 emu/cm^. 

At higher laser fluences with respect to those leading to Skyrmion generation (corresponding to an energy density 
up to 7 mJ/cm^ for each single laser pulse) we observed the formation of "doughnut" -shaped magnetic structures such 
as the ones shown in Figs. [2^ and [2)3. These structures are similar to those observed after cycling the magnetization 
of Co/Pt multilayers with perpendicular anisotropy [24]: the magnetization in the center is the same as in the 
surrounding medium, and is completely reversed in the annular region around the core. This texture can be viewed as 
two Skyrmions with opposite topological numbers (i.e. a Skyrmion and an antiSkyrmion) nested into each other. This 
results in a state with (5 = 0, corresponding to a nontopological soliton, which can be named Skyrmionium (Sk) in 
analogy with states of matter constituted by a particle bound to its antiparticle. Following this notation, a Skyrmion 
and an antiSkyrmion can be respectively abbreviated in Sk+ and Sk~ , since they can be obtained by "ionizing" Sk. 
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FIG. 2: (Color online) (a, b) Near- field Faraday rotation maps showing "doughnut" -shaped magnetic domains after single laser 
pulse irradiation (energy density ^ 7 mJ/cm^) in film areas showing opposite out-of plane magnetizations, (c, d) Corresponding 
Skyrmionium spin textures acccording to Eq. (|5]) (with Oo = 7r/2 for both Skyrmion and antiSkyrmion) . (e, f) Profile (dots) of 
the Faraday rotation measured along the lines in (a) and (b), respectively. The full lines in (e) and (f) are fits performed with 
the trial function describing Skyrmionium discussed in the text [Ro — 370 nm for the structure in (a) and Ro — 430 nm for the 
structure in (b)]. 



Sk is stabilized by the attractive force existing between a Sk+ and a Sk , generated by their opposite magnetic 
dipole moments. The simplest Sk spin texture with cylindrical symmetry is given by Eq. ([2j) with 



cr{r) = l-2e{r/Rofe-'^'^^o, 



(5) 



By substituting this trial function into our model Hamiltonian we are able to analytically estimate the value of the 
total energy Hx + Ha + i^D + Hz and to determine the Sk size Rq^ which is found to be 



Ro 



3.88(^ 



(4 - e)(7rr/^2 _ 2n/df) + AnAz/ea'^ 



(6) 



It is remarkable that, with no magnetic field (Az = 0), the ratio between the Sk radius Rq and the Sk+ radius Ri 
is a sample-independent quantity: Rq/Ri = 1.87. This ratio is in excellent agreement with the experimental value 
obtained from Figs. [1] and [2j Rq/Ri = 1.8 ±0.3. According to Eqs. (j4|) and (J6]), both Ri and Rq decrease with Q (the 
DDI strength), suggesting that a possible means to shrink the lateral size of such structures without the help of a 
magnetic field could be to employ such ferrimagnetic films at a temperature close to their magnetization compensation 
point, at which the two sublattices have opposite moments and the DDI vanishes. 

The evolution of the spin textures that are generated by single-laser-pulse illumination with increasing intensity 
is illustrated in Fig. [3l For fluences below 4 mJ/cm^ no magnetization reversal is induced, while single, isolated 
Skyrmions are generated in the 4-5 mJ/cm^ fluence range with 100% efficiency, i.e. each single laser pulse always 
gives rise to a Sk+. Conversely, fluences between 5 mJ/cm^ and 7 mJ/cm^ cause systematic Sk formation. The 
diameter of both Sk+ and Sk increases with fluence, suggesting that pinning defects might play a significant role in 
defining the ultimate size and shape of such structures. Apparently, Sk+ (Sk) is systematically generated when the 
laser fluence is over the minimum threshold [15| required to induce magnetic reversal over a wide enough area, with 
radius comparable to Ri (Rq). The role of defects might thus be to prevent the structures from collapsing to the size 
defined by Ri for Sk+ or by Rq for Sk. 

Upon illumination with laser fluences above 7 mJ/cm^, larger structures are generated, with topological properties 
that may vary from structure to structure. The texture shown in Fig. [3^ closely resembles to those observed in 
hexaferrite [9]. We interpret it as a combination of two Sk+ and one Sk~, resulting in a structure with Q = -\-l. It 
can be viewed as a di- Skyrmionium "cation" molecule Sk2", representing the topological analogue of di-positronium 
molecules. According to our model, concentric ring structures, made of alternating N (or A^ + 1) Skyrmions and A^ 
antiSkyrmions nested together, should all be stable. Their Skyrmion number is or ±1 depending on the number of 
Skyrmions being larger or equal to the number of antiSkyrmions. We find that the values of the zero-order parameter 
of adjacent ( ant i) Skyrmions should always be opposite, in agreement with the alternating helicity reversals inside 
Skyrmions with multiple-ring structure reported for hexaferrite [9;]. 
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FIG. 3: (Color online) Evolution of size and topology of magnetic domains generated in TbFeCo after single laser pulse 
illumination, as a function of laser fluence. The horizontal lines indicate the values of the Ri and Rq parameters evaluated 
for this TbFeCo sample (see text). The size is estimated by considering the rim of the domains defined by the condition of 
vanishing Faraday rotation and by averaging the radius of the domain over angle. The error bar on the size also considers the 
deviations from the axial symmetry of each considered structure. Insets: 2.5 x 2.5 /xm^ Faraday rotation maps. 

Figure [3F shows a different di-Skyrmionium texture consisting of one Sk+ surrounding two Sk~, yielding Q = —1. 
In such a Sk^ "anion" molecule the repulsion between the two Sk~ having parallel dipole magnetic moments is 
screened by the Sk+ that encircles them, which has an opposite dipole moment. It is still an open question whether 
Sk^ is intrinsically stable or a nonuniform distribution of spin-pinning defects in the polycrystalline film is essential 
to prevent this structure from collapsing or expanding. The non-centrosymmetric shape of the structures shown in 
Fig. [3] indeed suggests that defects are distorting the spin arrangement. Although we are not able to determine the 
conditions leading to the creation of either Sk2" or Sk^, their observation together with the reliable production of 
Skyrmionium represents an important step towards the study of complex topological states in condensed matter. 

Our results clearly demonstrate that by illuminating a thin TbFeCo film with single ultrafast laser pulses we can 
locally create stable magnetic configurations characterized by well-defined topological numbers. By tuning the laser 
fluence we can control the topology and create various bound states of Skyrmion molecules (Skyrmion "chemistry"). 
Such magnetic excitations are stable even without an external magnetic field. The ability of controlling the dimension 
and the number of Skyrmions or antiSkyrmions created during a single laser illumination event and the fact that such 
structures are topologically protected and thus intrinsically stable might have far reaching implications for the further 
development of non- volatile high-density magnetic recording and long term memories of the next generation. 
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